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Abstract 

We evaluate numerically the four-momentum-space gluon propagator in 
lattice Landau gauge, for pure SU (2) lattice gauge theory with periodic 
boundary conditions. Simulations are done in the strong-coupling regime, 
namely at (3 = 0, 0.8, 1.2 and 1.6, and for lattice sizes up to 30 4 , 16 4 , 20 4 and 
24 4 respectively. In the limit of large lattice volume, we observe in all cases 
a gluon propagator decreasing as the momentum goes to zero. This counter- 
intuitive result has a straightforward interpretation as resulting from the 
proximity of the so-called first Gribov horizon in infrared directions. 
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1. A non-perturbative investigation of QCD in the infrared limit is nec- 
essary in order to get an understanding of crucial phenomena such as color 
confinement and hadronization. In particular, the infrared behavior of the 
gluon propagator can be directly related []J to the behavior of the Wilson 
loop at large separations, and to the existence of an area law. Of course in 
developing non-perturbative techniques for QCD, one has to deal with the 
redundant gauge degrees of freedom. The standard gauge-fixing technique is 
based on the assumption that a gauge-fixing condition can be found which 
uniquely determines a gauge field on each gauge orbit. The correctness of 
this procedure was questioned by Gribov ||. He showed that the Coulomb 
and the Landau gauges do not fix the gauge fields uniquely, namely there 
exist gauge-related field configurations (Gribov copies) which all satisfy the 
gauge condition. Gribov copies affect also lattice numerical simulations [^]. 
In fact, the lattice Landau (or Coulomb) gauge is imposed by minimizing 
a functional which usually displays several relative minima (lattice Gribov 
copies). 

In order to get rid of the problem of spurious gauge copies, Gribov pro- 
posed @ the use of additional gauge conditions. In particular he restricted 
the physical configuration space to the region Q of transverse configurations 
(i.e. d-A = ), for which the Faddeev-Popov operator A4 [A] = — V- D [A] 
is nonnegative. This region is delimited by the so-called first Gribov hori- 
zon, defined as the set of configurations for which the smallest, non-trivial 
eigenvalue^ of the Faddeev-Popov operator is zero. We now know that Q is 
not free of Gribov copies and that the physical configuration space has to be 
identified with the so-called fundamental modular region 0, § . Nevertheless, 
the region Q is of interest in numerical simulations. In fact, configurations 
which satisfy the usual lattice Landau gauge condition belong to this region. 

The restriction of the path integral, which defines the partition function, 
to the region Q implies a rigorous inequality || |?J for the Fourier components 
of the gluon field. From this inequality, which is a consequence only of the 
positiveness of the Faddeev-Popov operator, it follows || 0] that the region 
Q is bounded by a certain ellipsoid G. This bound implies proximity of the 
first Gribov horizon in infrared directions, and consequent suppression of 
the low-momentum components of the gauge field, a result already noted by 

1 The Faddeev-Popov operator has a trivial null eigenvalue, corresponding to a constant 
eigenvector. 
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Gribov in reference . This bound also causes a strong suppression of the 
gluon propagator in the infrared limit. More precisely, Zwanziger proved 
that, in the infinite- volume limit, the gluon propagator is less singular than 
k~ 2 in the infrared limit and that, very likely, it does vanish as rapidly as 
k 2 . A gluon propagator that vanishes as k 2 in the infrared limit was also 
found (under certain hypotheses) by Gribov 0. In particular, he obtained 
D(k) = k 2 / (k A + 7). The momentum scale 7 appears when the restriction of 
the physical configuration space to the region Q is imposed. This propagator 
agrees with the zeroth-order perturbative prediction k~ 2 at large momenta, 
but gives a null propagator at k = 0. A propagator which is a generalization 
of the one obtained by Gribov has also been introduced in reference || as 
an Ansatz for a non-perturbative solution of the Schwinger-Dyson equation. 
Let us notice that these results are in complete contradiction with the k~ 4 
singularity obtained for the gluon propagator when the Schwinger-Dyson 
equation is approximately solved in the infrared limit [[J. 

The infrared behavior of the gluon propagator in lattice Landau gauge 



has been the subject of relatively few numerical studies [HJ [11], [T2|, [13|, [14 
In some cases |L2|, [14[ results seem to indicate that the gluon propagator is 
finite at zero momentum and in the infinite- volume limit. This is in agree- 
ment with Zwanziger's prediction of a gluon propagator less singular than 
k~ 2 . Recently, a gluon propagator evaluated avoiding Gribov copies |13[] has 
been successfully fitted by a Gribov-like formula jL5). Finally, in references 



lg, |17|, it was checked that the influence of Gribov copies on the gluon 



propagator (Gribov noise) is of the order of magnitude of the numerical ac- 
curacy. In particular, this seems to be the case even at small values of the 
coupling (3, namely in the strong-coupling regime, where the number of Gri- 
bov copies is higher and Gribov noise, if present, should be larger and more 
easily detectable. 



In references [16, 17 1 it was also observed, at (3 — 0, that the gluon 
propagator decreases as the momentum goes to zero. A similar behavior 
seemed to be present also at higher values of (3, and in particular at (3 = 1.6 
and for lattice volume V = 24 4 . In the last case, however, the behavior is not 
as clear as that at (3 = 0, due to the limited statistics (only 7 configurations). 
In the present work we extend the analysis done in references Jl6], [17| , and we 
would like to obtain a conclusive result for the infrared behavior of the gluon 
propagator in the strong-coupling regime. Let us notice that Zwanziger's 
predictions |5], [7| for the gluon propagator are (3 ' -independent in fact, they 
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are derived only from the positiveness of the Faddeev-Popov operator when 
the lattice Landau gauge is imposed. Thus, results in the strong-coupling 
regime are a valid test for those predictions. Of course, once the behavior of 
the gluon propagator is clarified in this regime, we should try to extend the 
results to larger values of /3, possibly up to the scaling region. 



2. We consider a standard Wilson action for SU(2) lattice gauge theory 
in 4 dimensions with periodic boundary conditions. For notation and details 
about numerical simulations we refer to [0, I7| . The only difference will be 



that here we consider lattices with different sizes in the different direc- 
Let us recall that the gluon propagator, which will be evaluated in 
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tions. 
lattice Landau gauge 



rt| |TSJ, is defined, in momentum space, 
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In Table [I] we report, for each pair (j3,N), the parameters used for the 
simulations. Only for the lattice volume V = 8 4 , at (3 = 0.8 and 1.6, and 
for V = 24 4 , at (3 = 1.6, we used the data reported in reference |L7|, average 
"fin"]. Computations were performed on the IBM SP2 at the Cornell Theory 
Center, on several IBM RS-6000/250-340 workstations at New York Univer- 
sity, on a IBM RS-6000/550E workstation at the University of Rome "La 
Sapienza" and on an ALPHAstation 255 at the ZiF-Center in Bielefeld. 

In Figure [l] we plot the data for the gluon propagator as a function of the 

2 Here k has components fc p = 0, 1, . . . , — 1. As in reference [[ItJ, we evaluate 
the gluon propagator by considering only values of k with three of the four components 
equal to zero, namely k = (0, 0, 0, k^). 

3 I thank S.Petrarca and B.Taglienti for kindly offering me the access to this machine. 
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square of the lattice momentum^ 



p 2 (k) =4^ sin 2 (vr k^) = 4 sin 2 (vr k 4 ) (3) 

for different lattice sizes, at (3 = 0,0.8,1.2 and 1.6 respectively.^ Finally, in 
Table |2|, we report the data for -D(O) at (3 = as a function of the lattice 
volume. 



Data at f3 = confirm previous results |L6], |17| : the propagator is decreas- 
ing (more or less monotonically) as k decreases. At j3 > the propagator is 
also decreasing in the infrared limit, at least for values of p 2 smaller than a 
turn-oveiQ value pf . Clearly p\ Q is /3-dependent. Using, for each (3, the re- 
sults for the largest lattice volume available we obtain: p 2 to w 1.4 at /3 = 0.8, 
p 2 to « 1 at = 1.2, and p 2 to « 0.6 at /3 = 1.6. 

As /3 increases the behavior of the gluon propagator depends strongly on 
the lattice volume. This is evident at (3 — 1.6. In fact, for the smallest lattice 
volume V = 8 4 , we obtain a gluon propagator increasing (monotonically) as 
k decreases. On the contrary, for V = 24 4 , the gluon propagator clearly 
decreases as p 2 goes to zero. 

Let us notice that, at high momenta, there are very small finite-size ef- 
fects, at all values of (3. The situation is completely different in the small- 
momenta sector, as already stressed for the case f3 = 1.6, and as can be 
observed in Figure |T| also for f3 = 0.8 and 1.2. The value -D(O) of the gluon 
propagator at zero momentum is also very volume-dependentQ (see Table |2] 
for /3 — 0, and Figure |I] for the other values of the coupling). In particular, 
at {3 = 0, as the volume increases, we obtain an initial monotonic decrease 
of the value of -D(O); for V = 18 3 x36 this value unexpectedly jumps up, and 
then starts again to decrease monotonically.^] Finally, for (3 = 0.8 and 1.6, 

4 As said in footnote we take k — (0, 0, 0, k^). 

5 For clarity, at (3 — we plot only the data for the three largest volumes; results are 
similar for smaller volumes. For the same reason, at /3 = 1.6, we do not show the data for 
the lattice volume V = 12 3 x24 which are very similar to the data for V = 16 4 . 

6 A drastic change in the behavior of D(k), around a turn-over momentum p to , has 
also been observed in reference |l4| for values of /3 in the scaling region. 

7 Only at = 1.2 the value of D(0) is essentially constant, within error bars, for the 
three lattice volumes considered. 

8 Notice that -D(O) gets the same value on different lattices (such as 16 3 x32 and 20 4 , 
or 20 3 x40 and 24 4 ) with "similar" lattice volume (respectively 131,072 and 160,000, or 
320,000 and 331,776). 
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the value of the zero-momentum gluon propagator decreases (more or less 
monotonically) as V increases. These results suggest a finite value for D(0) 
in the infinite- volume limit, in agreement with references \T2 , |T4] . From our 
data it is not clear if this value would be zero or a strictly positive constant. 
Therefore, the possibility of a zero value for -D(O) in the infinite- volume limit 
is not ruled out. 

We think that our data for the gluon propagator are very interesting. 
The prediction 0, ||, |7j of a propagator decreasing as the momentum goes 
to zero is clearly verified numerically for several values (3, even if only in the 
strong-coupling regime. Moreover, it appears that the lattice size at which 
this behavior for the gluon propagator starts to be observed increases with 
the coupling.^] Of course we should extend our simulations to higher values 
of f3 and to larger volumes. This is, at the moment, beyond the limits of our 
computational resources. 
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V 


config. 


therm, sweeps 


sweeps 


P 


0.0 


12 4 


100 


4 


2 


0.945 


0.0 


16 4 


100 


4 


2 


0.96 


0.0 


16 3 x32 


150 


2 


2 


0.97 


0.0 


20 4 


100 


4 


2 


0.97 


0.0 


18 3 x36 


150 


2 


2 


0.9675 


0.0 


20 3 x40 


150 


2 


2 


0.9725 


0.0 


24 4 


220 


4 


2 


0.975 


0.0 


28 4 


150 


2 


2 


0.97 


0.0 


24 3 x48 


90 


2 


2 


0.97 


0.0 


30 4 


70 


2 


2 


0.975 


0.8 


8 4 


200 


1000 


100 


0.875 


0.8 


12 4 


200 


1375 


125 


0.9 


0.8 


16 4 


96 


1100 


100 


0.94 


1.2 


12 4 


150 


1000 


100 


0.92 


1.2 


16 4 


100 


1000 


100 


0.93 


1.2 


20 4 


85 


1000 


100 


0.94 


1.6 


8 4 


200 


1100 


100 


0.865 


1.6 


12 3 x24 


200 


1650 


150 


0.91 


1.6 


16 4 


86 


1650 


150 


0.92 


1.6 


24 4 


34 


1650 


150 


0.92 



Table 1: The pairs (f3,V) used for the simulations, the number of configura- 
tions, the number of sweeps used for thermalization, the number of sweeps 
between two consecutive configurations used for collecting our data, and the 
parameter p used by the stochastic overrelaxation algorithm. 
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V 


D(0) 


12 4 


0.154(0.007) 


16 4 


0.149(0.006) 


16 3 x32 


0.143(0.005) 


20 4 


0.142(0.005) 


18 a x36 


0.160(0.005) 


20 a x40 


0.153(0.006) 


24 4 


0.153(0.004) 


28 4 


0.146(0.005) 


24 a x48 


0.144(0.007) 


30 4 


0.136(0.007) 



Table 2: The zero four- momentum gluon propagator D(0) [see eq. ([!])] as a 
function of the lattice volume V at (5 = 0. Notice that 16 3 x32 w 20 4 and 
that 20 3 x40 w 24 4 . 
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Gluon propagator; B = Q 



Gluon propagator; 8 = 0. 
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Gluon propagator; 6=1. 
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Gluon propagator; 6 = 1-6 
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Figure 1: Plot of the gluon propagator D{k) [see eqs. (|T|) and @j\ as a 
function of the square of the lattice momentum p 2 (k) [see eq. ([3])] for lattice 
volumes: (a) V = 28 4 (□), V = 24 3 x48 (O) and V = 30 4 (*), at (3 = 0; (b) 

V = 8 4 (□), K = 12 4 (O) and K = 16 4 (*), at j3 = 0.8; (c) V = 12 4 (□), 

V = 16 4 (O) and V = 20 4 (*), at /3 = 1.2; (d) 7 = 8 4 (□), V = 16 4 (O) and 
1/ = 24 4 (*), at (3 — 1.6. In all our runs we set k = (0, 0, 0, A^). Since we 
use periodic boundary conditions, only data for &4 < 1/2 are reported here. 
Error bars are one standard deviation. 
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